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Abstract 



The goal of this note is to construct a class of traveling solitary wave 
solutions for the compound Burgers-Korteweg-de Vries equation by means 
of a hyperbolic ansatz. A computational error in a previous work has been 
^ ■ clarified. 

00 

^ : 1 Introduction 

O ' Consider the following equations 

; o KdV: 

ut + auux + suxxx = 0, (1) 

o Burgers: 

X ' + "^^2; + I^Uxx = 0, (2) 

b : 

.•^ , o modified KdV: 

Ut + Pu^u^ + su^^^ = (s > 0) (3) 

These three equations play crucial roles in the history of wave equations. Equation 

(1) is named after its use by Burgers P for studying turbulence in 1939. Equation 

(2) was first suggested by Korteweg and de Vries [2j who used it as a nonlinear model 
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to study the change of form of long waves advancing in a rectangular channel. It was 
shown by Wadati |2i that equation (3) can be solved exactly by the inverse scattering 
method and the N-soliton solution can be expressed explicitly which is essentially 
the same as that in the KdV equation. A combination of these three equations leads 
to the so-called compound Burgers-Korteweg-de Vries (cBKDV) equation: 

ut + auux + l3u^Ur^ + ^Uxx + suxxx = 0, (4) 

where a, /9, ^ and s are real constants. The first term is the instationary term. The 
second and third ones are two non-linear convective terms with different orders. 
The fourth is the so-called viscous dissipative term: /i, which refers to a positive 
quantity, denotes the viscosity. The last one is the dispersive term. In order for the 
traveling wave to keep its appearance all through its propagation, an equilibrium is 
required between convective and dispersive terms. 

The coefficient a, which usually refers to the wave celerity, will be assumed positive 
in the following. The case a < can be easily deduced from the latter one by 
adimensionalization. 

The coefficient of the Uxxx term in the modified KDV equation being strictly posi- 
tive, we restrain our study to the case s > in (jlj). As in [8], the quantity —3(3 is 
supposed positive, which results in (3 < 0. 

In the past few years, the cBKDV equation has attracted continuous attention 
from a diverse group of researchers. Especially, the traveling wave solution has been 
explored extensively. An exact solution was presented by Wang using the ho- 
mogeneous balance method. Feng [5-7] investigated the same problem by applying 
the first integral method and the method of variation of parameters, respectively. 
A more general result was obtained by Parkes and Duffy jHl El using the automated 
tanh-function method and adapting the method of variation of parameters used in 
A generalization of equation (4) with higher-order nonlinearities in one- and 
two- dimensional spaces was treated by Zhang et al. (TUl HI] using the method of 
undetermined coefficients, and by Li et al. ^2] applying the tanh-function method 
and direct assumption method with symbolic computations, respectively. Describ- 
ing traveling waves of nonlinear evolution equations has been one of basic problems 
in theoretical and experimental physics. Traveling wave solutions to many nonlinear 
one-dimensional evolution equations can be derived from a set of ordinary differen- 
tial equations that can be interpreted as a flow in a phase space. In the last century, 
experiments on one-dimensional states of nonlinear traveling wave convection were 
undertook by many researchers in the narrow annular cell, condensed matter physics, 
plasma, particle physics, and so on [13-16]. 
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In the present work, we aim at extending results in a previous work [Sj, in which 
travehng wave solutions of the cBKDV equation were exhibited as combinations of 
bell-profile waves and kink-profile waves. Following [5], we assume that equation (4) 
has the traveling wave solution of the form 

u{x,t) = u{^), ^ = X — vt (5) 

where v is the wave velocity. Substituting(5) into equation (4), performing an inte- 
gration with respect to ^ and setting the integration constant to zero yields 

u"{i) + ru\i) + an(0' + &«(0' + cu{i) = 0, (6) 

where ^^=750 = ^, 6=^ and c = — ^. Equation (6) is a second-order nonlinear 
ordinary differential equation. We know that if all coefficients of equation (4) are 
nonzero, equation (6) does not pass the Painleve test [17, 18]. 

The rest of the paper is organized as follows. In Section 2, a class of traveling 
wave solutions to the cBKDV equation is presented by using a hyperbolic ansatz and 
the types of traveling waves are discussed theoretically. Numerical simulations for 
a couple of sets of parameters are illustrated. A calculational error in the previous 
literature is clarified at the end of this section. In Section 3, we give a brief discussion. 

2 Traveling Solitary Waves 
2.1 Hyperbolic Ansatz 

The phase plane analysis in Section 3 of Paper [5] provides us useful information 
when we construct traveling solitary wave solutions for equation (4). It clearly 
indicates that under certain parametric conditions, equation (4) does have bounded 
kink-profile traveling wave solutions. Based on this result, in this section, a class of 
traveling wave solutions of the cBKDV equation is established by using a hyperbolic 
ansatz which is actually a combination of bell-profile waves and kink-profile waves 
of the form 

n 

u{x, t) = ^ [Bi tanh* [d^x - vt + Xq)] + Di sech* [Ci{x - vt + Xq)]) + Bo, (7) 

i=l 

where the S-s, C-s D'^s, {i = I, ■ ■ ■ , n), v and Bq are constants to be determined, 
and Xq is arbitrary. 
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After substitution of (7) into equation (6), we get 



E 

i=l 



i (-sech[Ci + (! + «) sech[Ci 5,]'^+' sinti[zCi]'^)Cf Di 
+BiCf{{-l + i)i sech[Ci tanh[Ci 5]-^+' - 2i sech[Ci C]^ tanh[Ci ^f) 



II {-i sech[Ci £,]'^+^smh[Ci £,]CiDi+i Sech[Ci£,]'^BiCit!i.nh[zCi]''^+^) 

+— -s (8) 
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/3 (E (Si tanh'[C,fl+D, sech'[C, $])+Bo) 
J ml 

n n 

a (E (S> tanh'[C,5]+D, Sech'[C\ ^])+Bof v {B^ tanh"[C, sech'[C, ^D+Bo) 
+ 5 = 0. 

2s s 

The difficulty for solving this equation lies in finding the values of the constants Bi, 
Ci, Di, Bq and v by using the over-determined algebraic equations. Following jSj, 
after balancing the higher-order derivative term and the leading nonlinear term, we 
deduce n = 1. Then we replace sech(Ci^) by ^rrq^^^' sinh(Ci,^) by sElLi^Sll^ 

tanh(Ci ^) by and multiply both sides by (e*"^ ^ + e"*"^ ^) e^^ ^, so equation 

(8) can be rewritten in the following form: 

6 

J^PkiBo, B,, Ci, Di, v)e'''^^ = 0, (9) 

A;=0 

where the Pk {k = 0, 6), are polynomials of Bq, Bi, Ci, Di and v. 

With the aid of mathematical softwares such as Mathematica, when n = 1, equating 
the coefficient of each term e'^'^^^ [k = 0, ...,6) in equation (9) to zero, yields a 
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nonlinear algebraic system which contains seven equations as follows: 



Po 

Pi 

P2 

P3 

Pa 

P5 

Pe 



-2AvBo + 12aBo^ -8/3Bo^ + 24/iBiCi - 24 Di + 24 a Bo -Di +24/3Bo^-Di 
24 s Ci^ Di + 12 a + 24 /3 So i^i^ + 8 /3 -Di3 



-6i)Di + 6aBo Di +6/3Bo^-Di - 6aBi Di - 12/3 Bo -Bi Di 
+ e^Bi^Di +6/HC1D1 +6sCi2Di 



/ -9i;Bo + ^^i^ +3/3Bo3+3i;Bi -SaBoBi -3/3Bo2Bi - 

- a/JBoBi^ +3^Bi3 + 12/iBi Ci +24sBi Ci^ +6aDi2 
y + 12/3BoDi2-12/3BiDi2 

-12 1;^! + 12qBo -Di + 12 (3 Ba^ Di - 12 (3 Bi^ Di - 36sCi^ Di + 8 13 Di'^ = 0, 



0, 



- 9vBo 



'^"fo +3/3Bo^ -3^;Bi +3aBoBi +3/3Bo^Bi - 
3/3 Bo Bi^ -3i3Bi^ + 12/iBiCi -24sBi Ci^ +6aDi2 + 12 /3 Bq Di'^ 
12/3Bi Di^ 

-6i;Di +60B0D1 +6/3Bo^Di + 6a Bi Di + 12^Bo Bi Di 
6/3Bi^ Di -611C1D1 +6sCi'^ Di 



3^Bo + i2ii!_+/ji?(,3 
3/3BoBi2+/3Bi3 



3 Bi + 3 a Bo Bi + 3 /3 Bo^ Bi 



(10) 

System (10) can be solved consistently by using Mathematica again. 

For sake of simplicity, we use Si, £2, £3 and e to denote 1 or —1, and denote by n 

the quantity: 



a 



K — El 



The sets of solutions are given by: 



Bo 

Ci 



a. p 

20 ^3 



6 s 



£2 



6 s 



(11) 



£^3 

2k 

ie( 



(12) 



2 1/31 



S3 



6s 



'6s ^ 



6s ^* L2I/3I V 6s ^3 6sJ 4, 



with the constraint: 



or: 



£1 ^2 ^3 



(13) 



5i = Di = (14) 

Note that the second set to the fourth set of solutions are not of much interest, since 
they only correspond to three trivial cases — the constant solution. Now we examine 
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the first set of solutions. It is interesting to notice tliat tlie nature of traveling wave 
solutions completely depends on the values of the parameters a, s, /i. Moreover, 
we observe that 

i. depending on the sign of ^ + 2 s ~ ^^fsT ~ f^' traveling wave 

solution moves positively or negatively; 

a. the case of Si 7^ corresponds to a kink-profile wave solution; 

in. since —s /J > 0, Di is purely imaginary; equation (4) thus admits a complex 
traveling solitary wave solution. 

Due to the tanh function, the real part of the solution u has a kink-profile, 
while its imaginary part has a bell-profile (due to the sech function). 

2.2 Parametric Study 

In the following, we are going to show the influence of parameters a, /?, s, /i on the 
shape and velocity of traveling waves, and a couple of traveling waves corresponding 
to different sets of parameters are illustrated. 

2.2.1 Shape of Traveling Waves 

The variations of the coefficients Bi and Di have crucial influence on the bell-profile 
and kink-profile of traveling solitary waves. 



For this purpose, we consider the quotient of Bi and Di 



2. 



B,^ _ (2M-^i^2^7^) 

Due to the contraint (fT^ : 
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(15) 



It ensures: 



and, since Di is pure imaginary: 



£3 = 51 £2 (16) 

Bi^ = -Di^ (17) 



\Bi\ = \Di\ (18) 

which results in an equilibrium between the bell-profile and kink-profile of the soli- 
tary wave. 
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2.2.2 Velocity Analysis 



Here wc extend our attention to the wave velocity v. From the formula of v in the 
first set of solutions, we can see the rate of change of v (with negative option) with 
respect to a and n, respectively: 



da 



18\(3\ 



dv 
djji 



18 s\f3\ 

Since f3 < 0, there exists critical points for only one of those derivatives. For 
example, if £3 = —1, denote: 



V 6s 

In this specific case, we can derive variational tables as follows directly 



(19) 



a 


+00 


da 


- + 


V 


+00 

\ / 

4 s 





+00 


c)v 
dfj. 




V 


m 

\ 

—00 



Denote respectively etc the value of the parameter a which satisfies < ac, 
v{ac) = 0, and fic the value of the parameter ^ which satisfies < fic, v{fic) = 0. 
When a varies in [etc, +00 [ or /i varies in [0,^c]j the wave will propagate with a 
positive velocity. When a gets larger and larger alone, the wave will propagate with 
a big positive speed; but when fx becomes larger and larger alone, the wave will 
propagate with a big negative speed. 



Similarly, from the derivative of v with respect to /3 and s 
dv 



d\P\ 



36 /?2 



dv 



36 

s 



There exists critical points for only one of those derivatives. For example, if £3 = —1, 
denote: 



o? s 

6 



(20) 
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We can find the following table immediately 





+00 


m 


- + 


V 


\ / 

4s 



S 


+00 


Ov 

as 




V 


/ 

— 00 



As P approaches to zero, the absolute value of the wave speed will increase and 
eventually blows up. Similar thing will occur when s approaches to zero, but the 
wave will propagate with a positive speed. 



2.2.3 Numerical Example 

Figure ^ presents the real and imaginary parts of the traveling solitary wave as 
functions of the space variable x and the time variable t for a = 0.05, (3 = —0.15, 
s = 1, /i = 0.5. For this set of values: 



£1 — —£2 — — ^3 



(21) 




Figure 1: The real and imaginary parts of the traveling wave when a = 0.05, /? = —0.15, 
s = 1 and fi = 0.5. 

The real part of the traveling wave presents a kink-profile, while the imaginary part 
presents a bell-profile, as it could be expected, since |-Bi| = Di being purely 
imaginary. 

Here we wish to clarify that the coefficients of the solutions in [5, pp. 428-429] 
are incorrect. By comparison with the previous codes, we find that it is mainly due 
to Mathematica codes used for [5], in which a cubic term was missed while authors 
typed the algebraic system. 
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3 Discussion 



In this work, by applying a hyperbolic ansatz, we obtain a class of new kink-profile 
traveling solitary wave solutions to the cBKdV equation (1), which indicates the 
coefficients of the solutions in a previous work [5] are incorrect. One of the advan- 
tages of this approach is that it is very straightforward to generate new solutions 
and easily carried out with the aid of mathematical softwares. 

Although explicit forms of traveling solitary wave solutions are described herein and 
in the mentioned literature, as far as our knowledge goes, it is still unclear whether 
the cBKdV equation has any other type of bounded traveling solitary wave solu- 
tions. In the near future, we are going to use Lie group analysis as well as some 
other innovative methods to continue studying this challenging problem. Physical 
explanation for the solutions will be provided and the stability will be also analyzed 
as the parameters vary. 
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